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KOREWORD 


Satellite  altimetry  yields  geoid  height  data  along  the  ground  tracks.  These  data 
are  averaged  over  cells  into  which  the  earth's  surface  is  divided.  This  report  contains 
an  investigation  of: 

1.  The  error  committed  by  aggregation, 

2.  The  error  in  the  cell  average  due  to  sampling. 
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INTRODUCTION 


Satellite  altimetry  will  soon  yield  extensive  geoid  height  data  over  the  oceans. 
One  may  divide  the  surface  of  the  earth  into  approximately  square  cells  of  about 
equal  area  and  average  the  geoid  height  within  each  cell.  That  would  be  a method 
of  data  aggregation.  It  leads  to  an  imprecision  in  the  determination  of  the 
anomalous  field  from  geoid  height.  In  order  to  make  an  a priori  investigation  of 
this  effect,  one  may  make  the  following  hypothesis:  The  geoid  height,  without  the 
J,  term,  is  an  homogeneous  and  isotropic  random  field  over  a spherical  earth.1  - Its 
degree  variances  obey  Kaula's  rule.3  Random  quantities  are  underlined  in  this  report. 


THE  MATHEMATICAL  MODEL 


The  harmonic  expansion  of  the  geoid  height  is  written 


N(X.0)  = 


/.( m,C)U( m, t!:  A.0) 


where 


X is  latitude 
0 is  longitude 

The  functions  U(m,(d: \.<t>)  are  surface  harmonics  defined  by 


U( m, Id;  A.0) 


Cld+1) 


(Id  - lm|)I 
( C + |m|)! 


Plm  (sin  X )e,m 0 


for 


(d  = 0.1.  ••• 
m = -V.  -C+l.  — . C 
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These  harmonics  satisfy  the  following  orthonormality  relations. 

J"J"u(m',8’;  X.0H1E2  = 47r5(m.m':  8,8') 

They  have  been  defined  so  that  their  mean  square  over  the  sphere  is  unity. 

The  coefficients  in  the  expansion  of  the  geoid  height  are  random  variables 
satisfying  the  relations. 

I |Z(m,8) I = 0 


It)  1 0 R* 


( , , ) 10  " 
I <Z(m.8)Z<ni,8  )*  J ^ 


f>i  m, m';  8.8'  ) 


according  to  Kaula’s  rule  (R  is  the  nominal  radius  of  the  earth).  A discussion  of 
different  degree  variance  attenuation  laws  can  be  found  in  Rapp.  The  assumed 
homogeneity  and  isotropy  of  the  random  function  NtX.y'i)  defined  over  the  sphere, 
implies  that  the  harmonic  coefficients  Z(m,8)  are  uncorrelated  as  indicated. 

It  may  be  noted,  although  we  shall  make  no  use  of  it.  that  the  geoid  height 
autocorrelation  function  based  on  our  hypothesis  is  given  by 

, v oo  -ly  + ] 

I |N(X.<))N(X'.0')|  = It)  1 11  R:  I — 1\ ( cos  ) 


is  the  arc  lemith  between  the  rays  (X.0)  and  (X',0’). 


AVERAGING  OVER  CELLS 


The  earth's  surface  is  divided  into  zones  bounded  by  lines  of  constant  latitude 
at  equal  intervals.  The  zones  are  divided  into  equal  cells  along  lines  of  constant 
longitude.  A simple  consideration  shows  that  the  harmonic  expansion  of  the  error 
geoid  height  corresponding  to  these  cells  is 


w £ L(m,8)U(m,8;  X ,<p) 

d = n m = v’  ~ 


where 


L(m,t) 


v i zon'.njs  z _l!L 


«•=  2 111'  v' 


1 J A A(l) 


K(  l,J : m'A')*K<  I.J : mA) 


I and  J denote  the  zones  and  cells  within  a zone,  respectively.  A At  I ) is  the  solid 
ancle  subtended  bv  a cell  in  zone  1 at  the  centre  of  the  earth. 


The  quantities  Ktl.JimA)  are  defined  by 


M I.J ; tu.fc)  = — 

47T 


// 


U(m.V:X.0)*dS2 


coin  i.j  ) 


It  follows  that  the  error  anomalous  potential  is  given  by  (using  Bruns’  equation). 


AT  - 7 


l / R 


V + I 


1 = 0 m = - v \ r 


where  r is  the  radial  distance  and  y is  the  acceleration  ol  gravity  on  the 
earth. 


The  radial  error  gravity  disturbance  is 


6 Ac 


f)AT  y oc  , / |<V  + ; 
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£ ( V + 1)1 — ) L(mA’)U(ni.f;  X.0) 
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The  random  field  AN  defined  on  the  sphere  is  not  homogeneous  and  isotropic. 
Nevertheless,  in  view  of  the  orthogonality  of  the  harmonic  functions,  one  obtains  a 
simple  expression  for  its  mean  variance,  namely 


I j — ff  AN\lS2  | = S v pj|L(mA)l4 

I 47T  J J ) V*0  in  = - V ( J 


s p li  o r c 


The  expressions 
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I 
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o2<T;AN)  = 


« l s) 

£ EX  |L(m,C)|2  > 

ni=-li  f \ 


are  the  mean  degree  variances. 

Similarly,  one  obtains  for  the  mean  degree  variances  of  the  radial  error  gravity 
disturbance,  the  expression 


o2(C:H.A5gr)  = (B  + 1 )2 
Rw 


i 


j |L(m,8)| 


2j 
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where  r = R + H. 

Tlie  mean  degree  varianees  of  AN  are  a measure  - within  our  model  of  the 
violence  done  to  the  anomalous  field  by  the  operation  of  cell  averaging. 


NUMERICAL  RESULTS 


Calculations  have  been  made  for  a 5°  X 5°  partitioning  of  the  earth's  surface 
considering  harmonics  up  to  degree  60.  For  comparison,  the  corresponding  quantities 
pertaining  to  the  Kaula  field  have  also  been  calculated.  A recursive  algorithm  for  the 
integrals  of  the  associated  Legendre  functions  needed  for  the  calculation  of  the 
quantities  K(IJ:m.C),  was  devised  by  A.  R.  DiDonato.5 

In  Table  I we  present  the  degree  variances  and  accumulated  degree  variances  of 
the  error  geoid  height  and  the  nominal  Kaula  geoid  height.  The  factor  y 
(7  = 978  cmXs  2I  has  been  included  to  convert  these  into  potential  units.  The 
units  are  cgs. 

The  graphs  located  on  pages  I 2 through  44  of  this  report  present  the  logarithm 
of  the  error  potential  and  nominal  potential  at  the  surface  of  the  earth  cumulatively 
as  a function  of  harmonic  degree.  The  corresponding  logarithm  of  the  radial  gravity 
disturbance  at  different  altitudes  is  also  plotted. 


4 


One  observes  that  the  low  degree  variances,  of  .AN  are  small.  I lie  reason  is 
that  the  low  decree  harmonies  are  little  affected  by  the  5°  X 5°  cell  averaging,  l or 
higher  degrees  they  become  larger  and  finally  exceed  the  nominal  Kaula  values.  I he 
corresponding  harmonies  are  represented  with  an  error  exceeding  their  magnitude. 
However,  the  sum  of  the  degree  variances  of  AN  is  less  than  that  of  N ('•40-m  )• 
Asymptotically,  the  mean  degree  variances  of  AN  attenuate  like  const /ir  • I hose 
of  N like  const/n3. 


As  to  the  gravity  disturbance,  one  observes  that,  for  example,  at  an  altitude  of 
s()0  km  it  is  represented  to  a relative  precision  of  about  1,20.  At  low  altitudes  the 
low  degree  mean  variances  are  well  represented,  the  high  degree  ones,  poorly  as 
one  would  expect. 


DESCRIPTION  01  TABLE  1 


Column  1 The  degree  C 4 4 

Column  2 Mean  degree  variance  of  error  potential  in  units  of  cm  Is  . 1 he 

harmonics  are  so  normed  that  the  integral  of  their  square  over  the  unit 

sphere  is  47t. 

Column  3 Progressive  sum  ol  column  1 • 

Column  4 Degree  variance  according  to  Kaula's  attenuation  law. 

Column  5 Progressive  sum  of  column  4. 

7 = 478  cm/s2 
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Table  1.  Degree  Variance  of  5°  X 5°  Cell-Averaged 
Potential  and  of  the  Nominal  Kaula  Potential 


1 

tVi^an)  x io  9 

V 

y2o2(1t;N)X  10  9 

V 

0 

1 

. 1 99  X 1 0 

.175 

9703 

3 

. 1 93 

2875 

4 

.217 

1213 

5 

.243 

.828 

621 

14412 

6 

.270 

359 

7 

.296 

226 

8 

.323 

152 

9 

.349 

106 

10 

.374 

2.439 

77.6 

15333 

1 1 

.398 

58.3 

12 

.422 

44.9 

13 

.445 

35.3 

14 

.467 

28.3 

15 

.487 

4.658 

23.0 

1 5523 

16 

.507 

19.0 

17 

.5  26 

15.8 

18 

.544 

13.3 

19 

.561 

1 1.3 

20 

.576 

7.372 

9.70 

15591.8 

21 

.591 

8.38 

n 

.605 

7.29 

23 

.618 

6.38 

24 

.630 

5.62 

25 

.642 

10.457 

4.9  7 

15624.5 

26 

.652 

4.42 

27 

.662 

3.94 

28 

.671 

3.54 

29 

.681 

3.18 

30 

.689 

13.813 

2.87 

15642.41 

31 

.698 

14.512 

2.61 

32 

.706 

2.37 

33 

.715 

2.16 

34 

.723 

1.97 
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Table  1.  Degree  Variance  of  5°  X 5°  Cell -Averaged 
Potential  and  of  the  Nominal  Kaula  Potential  (Continued) 


1 7:u:(UN)  X 10  9 


35 

.732 

36 

.741 

37 

.752 

38 

.761 

39 

.773 

40 

.784 

41 

.800 

42 

.813 

43 

.831 

44 

.848 

45 

.871 

46 

.89 1 

47 

.9  1 9 

48 

.945 

49 

.980 

50 

1 .0 1 3 

51 

1.055 

52 

1 .097 

53 

1 . 1 50 

54 

1.203 

55 

1.271 

56 

1 .340 

57 

1.426 

58 

1.519 

59 

1 .633 

60 

1.764 

2)  -yVUi-.N)  X 10  9 

1 7.388  1 .8 1 

1 .66 
1.53 
1.41 
1.31 

: 1.149  1.21 

1.13 
1 .05 
.976 
.911 

25.361  .852 

.797 
.748 
.702 
.660 

30.109  .621 

.585 
.552 
.521 
.493 

35.886  .467 

442 
419 

398 
378 
359 


15653.33 


1 5660.47 


1 5665.38 


15668.907 


15671.525 


43.569 


15673.521 


DESCRIPTION  OF  THE  GRAPHS 


The  units  of  all  quantities  are  cgs.  except  that  height  above  the  earth  is 
reckoned  in  km.  The  degree  variances  correspond  to  harmonics  normed,  so  that  their 
mean  square  over  the  sphere  is  unity.  The  nominal  field  is  assumed  to  obe\  Kaula’s 
attenuation  law. 

Notation: 


N 

Geoid  height 

AN 

1 rror  in  geoid  height 

due  to  5°  X 5°  aggregation 

5 e 
- 1 

Radial  gravity  disturbance 

A6gr 

I rror  in  radial  gravity 

disturbance  due  to  aggregation 

II 

Height  above  Earth 

l..  f! 

Maximal  and  running 

degrees  of  harmonics 

o:(V  :AN) 

Mean  degree  variance 

of  AN  of  degree  t 

o:(t:N) 

Mean  degree  variance 

of  N of  degree  V 

a:(t:II.A6gr) 

Mean  degree  variance 

of  Afie  of  degree  V. 
- 1 

o’(?:ll,8gr ) 

Mean  degree  variance 

of  6e  of  degree  f 

|_r 

7 = 978  cm  s: 

I he  following  logarithmic  cumulative  degree  variances  are  plotted 

1.  log7‘w"(l  . ANl  = logK'  1 o’  (f:  AN I las  a function  of  degree  I. 

' 1=0  1 

2.  log  7"o’(  L:\ ) = log  { - «’(tl:N)|as  a function  of  degree  I 

1 i - n 1 


(ir.ipln  \l.irl  mi  \. 


X 


3.  log  o'(  L.  II:  A5gr)  = log  i X cr(l?:II.  A5g.)  J as  a function  of  degree  L and 
height  II. 


4.  log  o’<  L.  II;  6gf)  = log  { 1'  a2(t;ll,  6gr)  | as  a function  of  degree  L.  and 
height  II.  ‘v  = 0 ' ’ 


SAMPLING  ERRORS  DUE  TO  TRACK  SPACING 


Altimetry  yields  geoid  height  data  along  the  ground  tracks  of  the  vehicle.  We 
may  consider  the  following  simplified  problem.  A cell  on  the  earth,  bounded  by 
latitude  and  longitude  lines,  is  covered  by  a number  of  evenly  spaced  ground  tracks 
parallel  to  one  of  the  boundaries.  Along  these  tracks  geoid  height  is  given.  We 
estimate  the  mean  geoid  height  over  the  cell  by  averaging  these  data.  What  is  the 
variance  of  the  resulting  estimate  of  the  mean? 

Let  the  cell  be  defined  by 


-a  < 0 < a (longitude) 

-a  < A < a (latitude) 


and  let  there  be  k tracks  parallel  to  the  meridian  and  equally  spaced. 


The  mean  geoid  height  over  the  cell  is 


N 


_ V 

A v 


V 

III 


/.(m  A’) 


U(m,C;  X.0)  cos  XdXd 0 


where 


A is  the  solid  angle  subtended  by  the  cell  at  the  centre  of  the  earth.  The 
average  over  the  track  values  is 


A 

N 


I 

_ S' 

A 7 


S' 

111 


Z(m  A) 


U(mA:X  -cv  + (r-  ':)A)cosXdX 


One  obtains  after  a few  steps  the  following  expression  for  the  error  variance 


(N  \>- 


+ 1 


10  UIR- 


\- 


t (2t!)4  s o 


r I ( 2 1 + 1 I s n 


eu  f 
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(sin  X l eos  XdX 


\(  2s  + 1 ,2V  + I I I I’;  l (sin  X ) eos  XdX 


ftQ'.k 


f(a.k.2s)~ 

2s  + I ):  ! 


\(0A  > = \2V  + I 


N(niA)  = W2(2k  + h 


tv  md 
<k  + mi! 


ft  a . kan ) 


2sm  ma  2cv  sin  ma  ma 

— tor  not  a multiple  of  ir 

m k ma  p 1 

sin  -j- 

' -2a(-  I )rl  1 k>  for  m #0.  ™ = fir 


m 

0 


for  m = 0 


lere 


r is  am  integer 


One  in ;i\  easily  establish  that  I ■N:|  is  obtained  by  replacing  the  function  f by 


2 sin  ma 


c(a.k,m ) 


I 


m 


for  m i=  0 


for  m = 0 


It) 


I lie  results  calculation  s\  i 1 1 1 a 5X5  square  are  in  the  following  table 

A 

Number  of  Tracks  ()|N  N)  in  cm 


66 

18 

6.1 

2.7 

1.7 
1.14 

.84 

.41 

.10 


The  standard  deviation  of  N comes  out  to  be  38  m.  which  is  almost  the  full 
kaula  value,  because  of  the  relative  smallness  ol  the  cell.  We  see  that  a small 
number  of  well  placed  traeks  leads  to  a negligible  sampling  error.  The  principal 
sources  of  error  remain  the  uncertainty  in  the  satellite  ephemeris  and  the  sea  state. 


1 


4 

5 

6 
7 

10 

20 


Ri  l l Rl  NCI S 


1.  A.  \1.  Obukhov.  Statistically  Homogeneous  Random  lields  on  a Sphere.  U spell i 
Mat.  Nauk.  2.  No.  2 (1047).  pp.  166-168. 

2.  A.  M.  Vaglom.  Second-Order  Homogeneous  Random  lields.  Proe.  of  the  4'1' 
Bcrkclcv  Symposium  on  Mathematical  Statistics  and  Probability,  2.  pp.  563-622. 
Berkeley.  U.  of  California  Press  (1661). 

3.  W \1.  kaula.  Theory  of  Satellite  Geodesy.  Blaisdell  (1666). 

4.  Richard  II.  Rapp.  ( ieopotential  Coefficient  Behaviour  to  High  Degree  and  Geo  id 
Information  hy  Wavelength . Ohio  State  University.  Department  ol  (icodetie 
Science.  Rpt.  No.  180  (August  16^2). 

5.  A.  R DiDonato,  Recurrence  Relations  for  the  Indefinite  Integrals  of  the 
Xssociated  l.egendre  functions.  NSWC  l)L  leehnieal  Memo  No.  5 

( Publication  Pending). 


H = 5 Km 

LOG  <r2(L;H;  A 8*) 


I I i — I I I I I I I I r r T-r-r  r I I i i I i -1  1 I r 1 

21  24  27  30  33  36  39  42  45  48 

DEGREE  (L) 


T T T I » ' I 1 ' I ' 1 I r 
>1  24  27  30  33 

DEGREE  (L) 


' I ' ’ I ' 
36  39 


42 


’ I ' 

45 


H=  10  Km 

LOGo-2  (L,H,  A8gr) 


T rvT-i  i I i i I r i | i i \~r  i | i i | i ) '[ "i  r~T  i i I i i I 

27  30  33  36  39  42  45  48  51  54  57  60 

DEGREE  (L) 


LOG  <r2  (L,  H,8gr) 


■ M t 1 I r | I n T I I I I I r I 1 r i~t~t  i M t~T  i i 1 t -r-| 

27  30  33  36  39  42  45  48  51  54  57  60 

DEGREE  (L) 


18 


T ] T ~T— [ I I T T"l  | I I ( I I |~  I ~T  | 1 I [ T T | T "1  | I 1 | 1 T J T'TJ- | T ] 1 I | 
18  21  24  27  30  33  36  39  42  45  48  51  54  57  60 

DEGREE  (L) 


H = 250  Km 
LOG  0-2  (L,  H,  8gr) 
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H = 450  Km 

LOG  <72  (L,H,  A 8gr) 
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H = 600  Km 
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